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RESEARCH I^ORAKDUIg 

A PRSLUilHARY STUDY OF RAH- ACTUATED COOLING SYSTEMS 
FOR SUPERSONIC AIRCRAFT 
By Jacbson R. Stalder and Kenneth R. Uadleigh 


SUI£ I ARY 


An analysis has been made of the characteristics of 
several cooling cycles suitable for coclqpit cooling of super- 
sonic aircraft. All the cycles considered utilize the dif- 
ference between dynamic and ambient static pressure to 
actuate the cooling system and require no additional power 
source . 

The results of the study indicate that as flight speeds 
become greater, increasingly complex systems are required 
to reduce the ventilating air to tolerable temperatures. 

At altitudes above approximately 35,OCO feet, a system 
composed of an externally loaded expansion turbine In 
conjunction with a supersonic diffuser would maintain 
tolerable ventilating air temperatures, at least up to a 
flight llacli number of 2. The most complex system considered, 
composed of a compressor, intercooler, and expansion turbine 
with the Intercooler cooling air decreased in temperature 
by expansion through an auxiliary turbine is capable of 
maintaining a ventilating air temperature less than ambient 
temperature up to a flight Hach number of 3*7* The preceding 
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results for both systems are predicated on a cockpit pressure 
equal to ambient static pressure. 

It is possible that similar systems can be devised which 
will allow operation of ram-actuated cooling cycles with the 
cockpit pressurized, with, however, added system components 
required in the form of additional heat exchangers and turbines. 


INTRODUCTION 

It is generally realized that the problem of maintaining 
habitable cockpit temperatures in airplanes designed for 
supersonic flight will bo difficult, ; The necessity of 
cooling the pilot *s compartment has arisen in the operation 
of high-speed subsonic airplanes. The problem will naturally 
become much more acute as aircraft speeds are increasod 
through the transonic and into the supersonic speed range. 

The cooling problem in supersonic aircraft arises, in 
part, from the near-stagnation temperatures attained in the 
acceleration of ambient air to velocities approaching that of 
the airplane which prevail in the boundary layer. In effect, 
the airplane is surrounded by a thin layer of air at tempera- 
tures approaching stagnation value. Solar radiation into the 
cockpit through the canopy and the dissipation of heat by the 
pilot and by electrical apparatus further adds to the cooling 
load. In addition, the entering ventilating air which is at 
stagnation temperature must be reduced in temperature before 
admittance to the cockpit. 
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A refrigeration cycle utilizing atmospheric air as the 
working medium is currently being employed as a means of 
enclosure, cooling. This cycle uses air that has previously 
been compressed by the main engine compressor or by a cabin 
supercharger. The high-pressure high- temperature air is 
then cooled in an intercooler and expanded through a turbine 
to the enclosure pressure. 

At supersonic .flight speeds, the pressure rise occurring 
from the adiabatic acceleration of the ventilating air becomes 
of large enough magnitude that it may be possible to utilize 
the energy of the ram-compressed air to operate a refrigera- 
tion cycle as veil as for pressurizing the enclosure. 

It is the purpose of this report to examine the character- 
istics of several cooling cycles, which are actuated by the 
difference between dynamic or ram pressure and ambient static 
pressure, and to present the results of the analysis in as 
general a manner as possible. No predictions have been made 
concerning the magnitude of the enclosure cooling loads as 
this factor is dependent upon enclosure size, constructional 
details, amount of insulation employed, etc. Likewise, details 
of turbine and compressor speeds, sizes, types, have not 
been discussed in this preliminary report. 
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SY1IB0LS 

The following symbols are used throughout the report: 

Cp specific heat of air at constant pressure assumed 
constant at 0.2M- Btu per pound, °F 
e intercooler cooling effectiveness, dimensionless 

hp horsepower 

J mechanical equivalent of heat (77& ft-lb/Btu) 

M ilach number, dimensionless 

P stagnation pressure, pounds per square foot 

Q heat abstracted from air by intercooler, Btu per second 

T absolute stagnation temperature (°F + ^59*7) 

W weight flow rate of ventilating air, pounds per second 

W* weight flow rate of intercooler cooling air, pounds 

per second 

€ adiabatic shaft efficiency, dimensionless 

T[ duct efficiency, dimensionless 

*y ratio of specific heats of air (assumed constant at 1.40) 

Subscripts 


a,b arbitrary stations immediately upstream and downstream, 
respectively, from component under consideration 
c compressor 

1 ideal or theoretical 

o free stream 

t turbine 
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It* 



1 

1 * 

2 ; 

) stations as indicated in figure 1 
2 * ! 



ANALYSIS 


Five systems have been analyzed. Of these five, four 
are simple variations of the first basic system (System I) 

■which comprises a supersonic diffuser and an expansion turbine. 
In this basic system, the air is decelerated to zero velocity 
(relative to the airplane) in the diffuser and expanded from 
the resultant high pressure through a turbine to the pressure 
of the portion of the airplane being cooled — the enclosure. 

The turbine -warm is absorbed by an external load such as an 
electrical generator, hydraulic pump, or similar pioce of 
equipment , 

System'll is identical -with System I except for the 
addition of a heat exchanger between the diffuser and turbine. 
The heat exchanger employs an internal cooling medium, that is, 
fuel, liquid oxygen, or solid COg as the coolant. The turbine 
work is absorbed by an external load as in System I. 

System III employs the turbine work, to drive a compressor. 
The compressor is located downstream from the diffuser and 
increases the ventilating air pressure above the value of the 
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final diffuser air pressure. An air-to-air intercooler oools 
the air after compression and before it is expanded through the 
turbine to the enclosure pressure.’ 

System IV is Identical with. System III except that the 
Intercooler uses an internal cooling medium rather than free- 
stream air as a coolant. 

System V is similar to System III with, however, the 
addition of a second turbine which reduces the temperature of 
the intercooler cooling air. This secondary turbine Is directly 
connected with the main turbine-compressor system so that its 
work is delivered to the compressor. A schematic diagram of 
the five systems is sho T ,m in figure 1* 

For all systems analyzed, the simplifying assumption is 
made that the pressure drop of the air In passing through a 
heat exchanger is negligible in comparison with the pressure 
rise in the diffuser (and compressor when employed) and the 
pressure drop of the air in passing through the expansion 
turbine. The accuracy of this assumption is dependent upon 
the size and geometry of the heat exchanger and the validity 
of the asa\unption increases as the flight llach number Increases, 
Constant values of specific heats of air are also assumed. 

The air ‘Is’ assumed to be dry so that condensation effects are 

absent. ' ’ " ‘ 

The genoral relations employed in the analysis of the 
systems are as follows; 

The ideal isentropic pressure ratio resulting from 
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deceleration of air from an initial Ilach number 
zero velocity is given by 



“a 


to 


( 1 ) 


The adiabatic (not necessarily isentropic) tempera- 
ture ratio under the same conditions is 



7-1 
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M a 


3 


+ 1 


( 2 ) 


Equations (1) and (2) may be derived (reference 1) from 

consideration of the perfect gas law, the general energy 
equation, and the equation for sonic velocity in a perfect 


gas. 

A duct efficiency *n is defined as the ratio of the 
actual stagnation pressure rise obtained in the duct to the 
ideal isentropic stagnation pressure rise obtained for the 
same entrance llach number. Thus, 

_ (Pb — Pa) 

• (E*b “ ^a)i 

Rearranging equation ( 3 ) and substituting equation (1) 
there is obtained 



?b 

Pa 


1 + T| 
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2 




W 


The intercooler cooling effectiveness e is defined, in 
the usual manner, a& the ratio of the temperature drop of the 
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vcntilatlng air in passing through tho coolor to tho initial 
temperature difference "between the ventilating air and tho 
cooling air. 

The temporature ratio across a compressor or turbine is 
given, in terms of the pressure ratio, by equations ( 5 ) and 
( 6 ), respectively, 


« 1 + 


7-1 

(?b/P a ) 1 ~1 




In equation ( 5 ), e c is the adiabatic shaft efficiency 
defined as tho ratio of the isentropic temperature rise to the 
actual temperature rise of the air for. the. compressor pressure 
ratio. In equation ( 6 ), is similarly defined as the 

ratio of the actual drop in temperature experienced by the air 
as it drops in pressure passing through the turbine to the 
ideal isentropic temperature drop it would experience for 
the same turbine pressure ratio. 

The horsepower required to drive a compressor is given by 

h ' _ J=P w T a I" (Pb/Pa) 1 ^ - - 1 ] (7) 


and the horsepower delivered- by a ^turbine is 


J c p T a < t 
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The heat abstracted from the air "by an internally cooled 
heat exchanger is 


Q, = W c p (T a - Tb) (9) 

■which may be rearranged, in terms of the temperature ratio 


of the air in passing through the cooler to 


^b 

T q 


= 1 - 


Q 


^a 


( 10 ) 


The general scheme ’of analysis is to combine, suitably 
for each system, the preceding general .relationships in order 
to obtain the ratio of enclosure temperature to ambient 
static temperature in • terms of the flight liach number, the 
ratio of enclosure static pressure to free— stream static 
pressure and the efficiency of- the system components. 

A detailed analysis of each of the five systems is 

presented. in Appendix A* 

DISCUSSION 


For all systems, the final temperature ratio and the 
amount of enclosure pressurization that can be obtained is 
dependent upon the efficiency with which the diffuser converts 
the free— stream kinetic energy to static pressure. Avail- 
able quantitive information on the performance of supersonic 
diffusers is meager. In reference 2, Kantrowitz and 
Donaldson have prosented a method for the design of reversed 
DeLaval nozzlo-r-type diffusers and, in addition, have obtained 
test data over a limited range of liach numbers to check their 
analysis. The data of reference 2 was used in this report 
because of lack of comparable data on othor types of diffusers. 
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It is shown, in reference 2, that it is impossible to obtain 
a shock-free deceleration through the sonic velocity in a dif- 
fuser of this type and that an unavoidable loss in total 
pressure through the modium of a normal compression shock must 
result for the floitf to be stable. This loss inoroases with 
the entrance Ilach number. The data of Kantrowitz and Donaldson 
were used in this report for calculating diffuser efficiencies. 
The procedure used in determining the diffuser efficiencies 
used heroin was as follows: The maximum diffuser efficiency 

for a given entrance or flight 1-Iach number was calculated by 
the method of reference 2. This theoretical maximum efficiency 
was then multiplied by a factor, 0.93 > t 0 obtain an efficiency 
closely corresponding to the best test efficiencies obtained 
by Kantrowitz and Donaldson. It is worthy of note that, for 
a fixed geometry diffuser, the maximum efficiency occurs only 
at the design Ilach number. In this report, it is assumed 
that the optimum diffuser is used for each flight 1-Iach number. 
Flguro 2 shows the diffuser efficiencies used in this report 
as calculated by the preceding method. 

In order to visualize the magnitude of the stagnation 
temperature's which occur as a result of the acceleration of 
ventilating air, figure 3 has keen prepared. Figure 3 Is 
derived from equation (2), and shows the stagnation temperature 
as a function of altitude and Ilach number assuming HACA 
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standard air properties. (See reference 3») 

The maximum degree Of enclosure pressurization obtain- 
able from flight ram, shown in figure 4 -, was calculated by 
means of equation ( 4 ) using diffuser efficiencies taken from 
figure 2, Constant pressure lines for several effective 
enclosure altitudes are also shown in figure 4 , The lower 
limit of the curves was taken at 3 0,000 feet altitude, since 
it is unlikely that prolonged supersonic flight would be 
undertaken below this altitude. The cooling systems discussed 
herein would operate less effectively as the enclosure 
pressure increased and would not operate if the enclosure 
were maintained at the maximum possible ram pressure. It 
is probable, however, that systems of this type could be 
devised which would allow almost complete ram pressurization 
to be utilized and still maintain the system effectiveness, 
with, however, the addition of more pieces of equipment - 
turbines, heat exchangers, etc. 

It is unfortunate that the equations for the temperature 
ratio across each system do not lend themselves to plotting 
in terms of nondimensional or dimensional groups of variables, 
so that the effect of a change in efficiency of a system 
component is immediately apparent. In order to calculate 

the performance of the system it is necessary to assume 
values for each of the component efficiencies. The follow- 
ing numericaT values wore used for substitution in the equa- 
tion for the tomperatwe ratio across each systemt 
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1* Diffuser 'efficiency, taken from figure 2 

2. Turbine adiabatic efficiency, O.S 

3. Compressor adiabatic efficiency, 0.7 
4-. Air-to-air intercooler cooling effectiveness, 0.9 

The enclosure ws,s assumed to be unpressurized, hence the 
enclosure pressure was taken as equal to ambient static pressure. 

It is thought that the above values approximate the 
maximum efficiencies that are practically obtainable, consider- 
ing the probable small size of the equipment. The seemingly 
high value of intercooler effectiveness arises from the fact 
thab the high ram pressures allow the use of multipass (and 
hence high effectiveness) heat exchangers. 

The performance of System I is shown in figure 5 . The 
pertinent point concerning this sytem is that the final 
temperature ratio is always greater than unity, that is, the 
entering ventilating air temperature Is always higher than 
free-stream static temperature, owing to energy losses in the 
diffuser and turbine. 

The temperature ratio of the air after passage through 
System II is shown in figure 6. In figure 6, the parameter 
Q/WctjTo represents. the fraction of the initial total heat 
content of the ambient air that is removed by the Intornal 
oooler. This variable is a function of the effectiveness of 
the internally cooled heat exchanger as well as the lowest 
temperature obtainable from the cooling medium",' From figure 6, 
it is apparent that if a suitable cooling medium can be 
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employed. It is possible to obtain very low outlet tempera- 
tures from this system. 

The determination of the oharact eristics of System III 
involves , as noted in the appendix, the graphical solution of 
equations (A9) and (A10) in order to determine the compressor 
pressure ratio values for subsequent substitution into equa- 
tion (AS). The values of the compressor pressure ratio 
determined in the foregoing manner are shown in figure 7« 

The performance of the system as represented by the final 
temperature ratio of the ventilating air is presented in 
figure S. It is noted that this system is capable of main- 
taining a system outlet temperature loss than free-stream 
static temperature up to a flight Ilach number of approximately 
1.7. 

Tho performance of System IV Is dependent upon the amount 
of heat abstracted from the air by the internal cooler, as in 
System II. The compressor pressure ratios obtained from a 
graphical solution of equations (A13) and ( Al4) are shown in 
figure 9» She pressure ratios decrease rapidly as heat Is 
abstracted by the Internal cooler due to tho decreased 
turbine worh available with the lower turbine Inlet tempera- 
tures. The final over— all temperature ratio across the system 
is shown In figure 10 . It may bo soon from this figure that 
the final temperature ratio Is determined by the amount of 
heat abstracted by the internal cooler and that tho system 
is sensitive to changes in flight 1-Iaoh number since a slight 
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change in llach number increases the final temperature ratio 
considerate . 

The determination of the performance of System V is 
complicated by the presence of the auxiliary cooling turbino 
that decreases the temperature of the cooling air before its 
passage through the intercooler. It becomes necessary to 
assume values for the mass flow ratio of intercooler cooling 
air to ventilating air, as may be seen from an examination of 
equation (A1S ) . In order to obtain the system temperature 
ratio as repi-esented by equation (A17) it was necessary to 
solve, graphically, equations (AlS) and (A19) for values of 
the compressor pressure ratio. The numerical, values obtained 
are shown in figure 11. From a consideration of figures 11 and 
12, it is evident that' the compressor pressure ratios become 
extremely high at llach numbers greater than 2.5 for a value 
of h' 1 /V7 of 2. It will be noted that these pressure ratios 
are considerably hi^ier than those obtained in Systems III 
and IY, due to the added work put intro the compressor by the 
auxiliary turbine. The auxiliary turbine performs two 
beneficial functions — it increases the work available to 
drive tho compressor' and, henco, inoreasos the pressure ratio 
across the main turbine, and it provides cold exhaust air to 
cool tho primary ventilating air. 

The final vent i la ting-air- temperature ratio across System 
Y is shown in figure 12. The performance of this system shows 
a marked decrease in over— all temperature ratio compared to 
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that of the other systems due to the "beneficial effect of 
the auxiliary turbine. From figure 12 it is seen that, for 
a value of mass— flow ratio of cooling to ventilating air of 
1.5, the system outlet temperature may bo maintained less 
than, or equal to, ambient static- temperature up to ,a llach 
number' of 3.7* The added performance Of this system is, of 
course, accompanied by an increase in tho power required 
' to ram the additional cooling air through the system. 

The fact that the temperature ratio decreases to a 
■ minimum- at a llach number of approximately* 1.9- and .then 
increases' as the ‘flight llach number increases is -due to the 
■opposing effects of increasing pressure ratio across, the aux-* 
•iliary" turbine- and •‘decreasing duct efficiency,.-. 

For all the systems ■ discussed, it is strongly . empha- 
sized that the performance shorn is the maximum obtain- 
able with. ‘-the. assumed values of compressor and.. turbine 
efficiencies and intercooler effectiveness, because the 
maximum obtainable value of diffuser efficiency for the L type 
of diffuser considered was tamcn for each flight 1-Iach number. 
The performance over a range of " fligjit llach numbers of any 
"actual enclosure cooling system with a fixed geometry diffuser 
of the type discussed would decrease below the ideal values 

I 

shown heroin. The values of diffuser efficiency used herein 
should not be considered the maximum obtainable with any type 
of diffuser. It is quite possible that diffusers can be 
designed mien will have performances superior to the simple 
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reversed DeLaval— nozzle type discussed. 

Systems, such as those described in this report, which 
take free stream air and pass it through a cooling cycle con- 
tribute drag to the airplane by virtue of the difference in 
momentum between the air taken into, and that discharged from, 
the airplane. Thus, the choice of a particular system will 
depend upon the system drag characteristics as well as upon 

the internal efficiency" of the installation. The drag 

* 

produced by a system is greatly dependent upon the size end 
configuration of the ducts and heat exchangers, hence it is 
not possible to draw specific conclusions concerning the drag 
of the several systems discussed herein. Hovever, it is 
apparent that cooling systems that handle relatively large 
amounts of cooling air such as Systems III and Y will have 
unfavorable drag characteristics and it may be necessary to 
reduce the heat-exchanger effectiveness in order to eliminate 
excessive pressure drops on the cooling air side. Detailed 
computations required to design an optimum system would require, 
as a starting point, specifications concerning flight speeds, 
altitudes, and cooling load, and hence have not been undertaken 
in this general report. 

CONCLUDING- RSIIAHKS 

The high ventilating and boundary-layer-air temperatures 
attained at supersonic velocities will probably require use of 
an enclosure cooling system to permit piloted supersonic aircraft 
to be flown. 


NACA Eli No. A7G04- 


If 


An examination of the characteristics of several ram- 

actuated cooling ^systems discloses that these systems may, 

from a standpoint of cooling performance, be suitable for the 

cooling of aircraft or missiles operating at supersonic 

velocities. The operation of all the systems depends upon 

the difference in pressure • between the enclosure or cockpit 

and the ram pressures due to the velocity of the aircraft. 

A system (System V) composed of a supersonic diffuser, a 

compressor, intercooler, and two expansion turbines appears 

promising, from a cooling standpoint, in view of the fact that 

the ideal performance of the system indicates a system 

outlet temperature less than ambient static temperature up 

to a flight iiach number of 3*7; provided the enclosure is 

maintained at ambient static pressure. 
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APPENDIX 
System I 


For this basic systora, equations (2), (4), and (6) are 
combined to give the final temperature ratio across the 
system in terms of the initial (flight) llaoh numbor, diffuser 
efficiency, turbine efficiency, and the ratio of enclosure 
static pressure to frec-stream static pressure. 

The latter ratio is, for a given system and flight 
liach number, dependent ' upon the amount of enclosure pres* 
surization desired. Combination of equations (2), (4), and 
(6) gives, for the final temperature ratio across the system. 


^ = . (Al) 

To VTo/V^/ 









(A2) 


The powor available from the turbino, which is absorbed by an external load, 
is, 
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System IX 

The over-all final temporaturo ratio across the Systran Is givon by, 


(A3) 
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Substituting values for and there is obtained 
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It will bo noted that equation (AS) contains tho term Ps/Pi, the compressor 
pressure ratio* This torn may bo ellmlnatod by equating the expressions for compressor 
and turblno powers since all of tho turbine powor is delivered to tho compressor. 

Honoo, whon equations (7) and (S) are oquatod and the rosultant expression Is solved 
for the term P a /Pi tho following quadratic solution rosults: 



Values of tho term T1/T3 in equation (A 9 ), are obtained from the follow- 
ing expression which is derived from oquation (5) and the expression for intor- 
coolor effectiveness, 
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To obtain values of the term P 2 /Pj, for substitution into equation (A3), 
a graphical simultaneous solution of equations (A9) and (A10) Is made* Values 
of the term P^/Pq may be chosen for any degree of enclosure pressurization 
desired* 

System XV 

13m} final temperature ratio across the. system is, 


T* 

T « 


T 0 ATr 

or. In terms of tho system variables, 


^ To AiAt-At./ 
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The unknown compressor pressure ratio P a /P x , in equation (A1H), is determined 
in the some manner as for Systom III. The resulting equations involving Pg/Pj, are 


ro 

v>l 
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Equations (JU/j) and (£L^) arc simultaneously solved for values of P3/P1 
whioh are subsequently substituted into equation (il2). It will be noted that 
in this sytem, an additional independent variable Q/WopT 0 is introduced* This 
variable is determined by the amount of heat abstracted from the air by tho inter- 
coolor and may bo arbitrarily olios on to yield the dcslrod final temperature ratio. 


System V 
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Tho analysis of this systom is similar to that of System III except that 
the intorooolor oquatlon must bo modified as the cooling air is now precooled 
by expansion through a turbine. The expression for intercooler effectiveness 
is, for this systom, 


= T 3 - T 3 
Ta - T 3 1 


(A15) 


The assumption is made that the preoooling turblno has the same adiabatio 
efficiency as the main turbine. The expression for the over-all final tempera- 
ture ratio across the system then is 
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In terms of the system variables, 
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(A17) 

where the term P^Pq is taken from equation (4-), In the derivation of equation 
(A17) it is assumed that P 3 l is equal to Po, that is, the intercooler pressure 
drop is considered negligible. 

The unknown - compressor pressure ratio Pa/P* in equation (A17)> determined, 
by equating the compressor power to the sum of the powers of the main and precooling 
turbines. The resulting expression for the term Pa/Pj la 
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In the simultaneous graphical solution of equations (2S) 
and (29), it is necessary to assume values of the mass-flow 
ratio of intercoolor cooling air to ventilating air W /W. 

In the actual case, the value of this parameter is detexmilned 
by the area of the intercooler heat-transfer surface and the 
over— all heat-transfer coefficient obtained in the cooler, 
since the cooling effectiveness is arbitrarily chosen. 

The presentation of the results of the foregoing analysis 
in the form of a final system temperature ratio is believed to 
be the most convenient form of presentation because tempera- 
ture changes due to variation in altitude then do not enter 
the equations . 

The cooling capacity of any system may bo calculated for 
any desired flight Ilach number and altitude by the following 
method : 

The ambient temperature is multiplied by the system 
temperature ratio to obtain the temperature of the ventilating 
air entering the enclosure. The difference between the enclo- 
sure ambient temperature and the entering ventilating air 
temperature is then, the temperature difference available for 
cooling tho enclosure, llultiplication of the temperature 
difference by tho term Wc p gives the oooling capacity of 
the system in terms of Btu per second. 
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